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Consider the equation 


x = X(x). (1) 

Assume that the behavior of the family of solutions is known, expressed by the estimate 

lx(t,x^)-x(t,x,) I < Kq f(t) lx, -Xg I , t >0, 

Here K„ is constant, and f is normalized by f(0) = 1 , f (t) > 1 . The function f (t) then 
expresses the stability or instability of ( 1 ), so that stability is equivalent to f = 1 , and a 
typical instability would be f (t) = t + I . 

Suppose we wish to numerically compute a particular solution of (1), say x (t, x^), for some 
fixed x„ . A remarkable theorem of Babuska asserts (under reasonable assumptions on X) 
that if a strongly stable difference scheme is used, numerical integration of (1) with x (0) = x„ 
is equivalent to obtaining an exact (theoretical) solution of 

y = X(y) + G(t,y), (21 

where y (0) = x, , (here lx, - x„ I is the starting error), and I G (t, y) I < rj, (here r) (small) 
reflects the discretization and truncation error). 

Babuska used this result to obtain global error bounds for (2), assuming x (t, x„) is a 
uniformly asymptotically stable soluvion of (1). Here we wish to compare solutions of (1) 
with solutions of (2), with particular attention given to the effect of different types of 
f(t). 

To develop the comparison, a converse theorem using Lyapunov functions due to Yoshizawa 
and Hale is applied. Let z = x - x (t, x^ ), then ( 1 ) becomes 

z = Z (t, z), (3) 

and the new estimate becomes 

• z (t, t^, z„) - z (t, t„, z,) I < K„ f (t) I z, - Zo' 
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f(0 , 

= K (tg) I z, - z^ I exp (a (t) - a (t^)) 


( 4 ) 


where K (t^) = f (t^), a (t) = log f (t). 

Then one can show there exists a function V (t, z) satisfying 

I zl < V (t, z) < K (t) I zl . 

I V (t, z) - V (t, y)l < K (t) I z - yl , and 


(a) 

(b) 


V 




(c) 


where Vj denotes the derivative of V along solutions of (3). Note that (a) and (c) allow 
one to recover (4), using standard theorems. Further, it is easy to show, using (b) that 


Vj <iv + K(t)n= Y 

Then an immediate consequence is that 

I y (t, X,) - X (t, x„)l < K„ f (0 1^1 X, - x„l + T?tj 

So we see that the term Kj, f (t) I Xj - I reflects the propagation of the starting error, 
and the term r? f (t) gives an estimate of the effects of the discretization and truncation 
error. 

Evidently, it is better to integrate a stable system, f = 1, than an unstable system, 
f (t) = t + 1 . A simple example is the unperturbed two-body problem, which is unstable, 
as the difference between two periodic solutions of different periods grows (at least locally) 
like t + 1. An example of a stable problem is the stabilized Kepler problem, as given by 
Baumgarte and Stiefel, where now all solutions have the same period. 

To extend the above to perturbed problems, consider 

i = X (x) + eX, (t. X) (5) 

Assume that in some region R, I X, (t, x) I • ' M, Here one supposes this region R contains 
the solution of interest and is large enough so that continued integration would not be of 
interest long before the integrated solution leaves R. 
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This assumption will allow the effect of f (t) on the error to be observed, for now we 
consider 


y = X (y) + eX, (t, y) + G (t, y). (6) 

and y (0) = X, , I G (t, y) K t? as before. 

The same V-function leads to the estimate 

I y (t, X,) - X (t, Xj,)l < K„ f (t) 1 ^ I X, - x„l + t (tj t eM)j , 

which is valid as long as both solutions remain in R. Ar^ain one can see the different error 
estimates obtained for f = 1 o*" f = t + 1 . Remark ; The foregoing gives upper bounds on 
the error. The strength of the conclusion then rests on the sharpness of the upper bounds, 
for one cannot conclude Error, < Error^ on the basis that a crude upper bound of the first 
is less than a crude upper bound of the second. 

However, by considering simple examples, one can see that at least in these cases, numerical 
integration of a stable system gives an error growth 0 (rjt), while for an unstable problem, 
the error grows like 0 (rjt^ ). So to that extent the foregoing estimates seem reasonable. 

Further, numerical experiments are being conducted at Goddard Space Flight Center using 
the stabilization techniques of Baumgarte and Stiefel, and the results should give further 
insight into the significance of the above results. 

In this regard, perhaps it should be observed that a similar approach to obtain lower bounds 
does not seem possible, as the difference cf two almost periodic functions will, in general, 
have an arbitrarily small lower bound. 

Note that the above . .imates can also be obtained in more complicated problems, with 
such means as time-dependent potentials. Details will appear in a forthcoming paper. 
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Optimal explicit Runge-Kutta methods are developed for solving the initial value problem 
for systems of ordinary differential equations. These methods have an optimal estimate of 
the local truncation error term, thereby allowing the option of implementing a variable- 
step strategy. The coefficients of the Rui^e^Kutta method are selected so that the local 
truncation error is minimized and that the absolute stability region is maximized. 
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